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The  magnitudes  of  the  experimental  and  calculated  responses  for  the  three  conditions 
(V,  =  0.  V,  =  0,  V,  =  V,)  are  shown  in  figures  5,  6,  and  7.  As  expected,  the  agreement  right 
at  the  resonance  points  is  poor  because  of  the  lossless  line  assumption.  However,  the  comparisons 
of  experimental  and  calculated  results  over  the  remainder  of  the  curves  seem  to  he  more  than 
sufficient  to  validate  the  applicability  of  superposition. 
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I.  INTRODUCTION 

In  analyzing  the  EMP  response  of  complex  systems,  the  analyst  is  usually  forced  to  make 
simplifying  assumptions  about  the  system  and  its  electrical  characteristics.  If  the  interaction  and 
coupling  problem  of  interest  is  an  exterior  cable  (one  exposed  directly  to  the  incident  EMP)  it  is 
normally  necessary  to  assume  that  no  other  cables  associated  with  the  system  are  excited  by  the 
EMP,  that  the  cable  has  a  constant  height  above  ground  (or  below),  and  that  the  cable  is  straight 
with  terminations  appearing  only  at  each  end.  In  many  instances,  these  necessary  assumptions 
are  not  unduly  restrictive  and  some  worst-case  “hand- waving”  can  be  done  to  justify  making 
them.  However,  these  justifications  become  extremely  difficult  for  complex  systems  which  have 
many  cables  associated  with  them. 


This  paper  investigates  a  technique  for  extending  the  applicability  of  existing  single-cable 
coupling  codes  to  complex  networks  of  uncoupled  cables.  The  technique  uses  the  property  of 
superposition  to  account  for  the  effects  of  having  an  arbitrary  group  of  cables  connected  together. 
The  solution  presented  here  is  limited  to  the  response  of  a  cable  which  has  an  arbitrary  number 
of  other  cables  attached  to  its  ends  as  shown  in  figure  I.  No  loops  of  cables  are  permitted.  The 


A  system  of  m  connected,  uncoupled  cables  excited  by  an  electric  field  E  propagating 
in  the  Y  direction. 


extension  of  this  technique  to  include  loops  would  require  a  two-port  distributed  system  repre¬ 
sentation  of  all  the  cables  involved  in  the  loop;  this  representation  does  not  appear  to  have  a 
compact  solution  at  this  time.  The  present  technique  is  being  used  to  modify  the  transmission¬ 
line  coupling-code  program  FREFLD  and  the  results  of  this  effort  will  be  reported  on  in  a  future 
report. 

2.  APPROACH 

If  one  conducted  a  review  of  the  various  frequency-domain  solutions  for  the  current  response 
of  a  two-wire  transmission  line  excited  by  a  horizontally  polarized  plane  wave,  most  of  these 
solutions  would  be  similar  to  the  equation1 


v  £(»)  f 
,U)--20  { 


[A  cosh  f  (/  -  x)  +  Zt  sinh  f  (/  -  *)] 


f(K  -  Z0)[l  ~  *-«•♦*’*]  (K  +  Z0)[ I  - 

x[  v  +  jp'  r  -  jp1  ) 

+  e~ie',(K  cosh  rx  +  Z„  sinh  Tx) 

/(K  -  Z,)[l  -  *«r-*xr-i>]  +  Z()[i  _  ■ 

x[  r  -jp'  r  +  jp' 


j(K-Z,)[  1  -  e'r 

\  r  -  jp- 


where 


£(o>)  is  a  source  term  relating  the  incident  field  to  a  voltage  exciting  the  transmission 
line, 

D  =  A[(A*  +  Z„Z,)  sinh  VI  +  K(Z0  +  Z,)  cosh  17], 

K  is  the  impedance  of  the  line, 

T  =  a  +  jp  is  the  propagation  characteristic  of  the  line, 
a  is  the  attenuation  constant  of  the  line, 

P  =  oi/V  is  the  phase  constant  of  the  line, 

V  is  the  propagation  velocity  of  the  line, 

I  is  the  line  length, 

P '  =  at  sin  ( 0 )  cos  (<li)IC  is  the  delay  associated  with  the  plane  wave  propagating  from 
one  end  of  the  line  to  the  other, 

<o  =  2ir/(where  /is  frequency), 

0  =  angle  between  the  electric  field  and  the  axis  of  the  cable, 

ifi  -  angle  of  incidence  between  the  field  propagation  vector  and  a  plane  normal  to 
the  line, 


’  R.  F.  Gray.  Nuclear  Electromagnetic  Pulse  Simulation  hy  Point-Source  Injection  Techniques  for  Shielded  and 
Unshielded  Penetrations.  Harry  Diamond  Laboratories.  HDL-TR-I7J7  (December  I V75). 


c  =  speed  of  light, 

x  =  point  at  which  the  current  is  calculated, 
Z0  =  termination  at  x  =  0,  and 
Z(  =  termination  at  x  =  /. 


The  geometry  used  is  shown  in  figure  2.  This  solution  was  developed  for  a  wire  or  cable 
over  a  ground  plane  as  shown,  but  the  geometry  only  affects  the  source  term  E(a>)  and  not  the 
rest  of  the  equation.  The  solution  for  excitation  by  a  vertically  polarized  plane  wave  has  a  form 
similar  to  that  given  here;  the  only  differences  are  in  the  trigonometric  terms  and  the  grouping  of 
terms  in  the  numerator.  Therefore,  it  will  be  stated,  but  not  proved,  that  the  technique  developed 
here  can  be  equally  applicable  for  vertical  and  horizontal  polarizations. 
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In  order  to  solve  the  coupling  problem  shown  in  figure  I,  we  first  assume  that  each  of  the 
m  -  I  cables  connecting  to  cable  n  can  be  represented  by  its  own  equivalent  source, 
which  is  either  a  current  or  a  voltage  source,  as  shown  in  figure  3.  Now,  if  the  equivalent  sources, 
5*.  exist  and  the  only  interaction  between  the  cables  is  at  x  =  0  or  x  =  /  but  not  both,  then  it 
should  be  possible  to  develop  a  systematic  approach  to  solve  for  the  current  of  interest,  /„(*), 
at  any  point  x,  using  the  principal  of  superposition.  It  may  be  possible  to  account  for  the 
interaction  between  the  cables  by  introducing  additional  equivalent  source  terms.  There  would 
also  be  an  effect  on  the  cable  impedance  if  there  is  an  interaction  between  cables  over  their 
length.  In  this  report,  it  is  always  assumed  that  the  areas  of  interaction  between  cables  are  small 
compared  with  the  areas  of  noninteraction.  This  limitation  does  not  appear  to  be  overly  restrictive, 
since  a  great  many  configurations  are  not  in  conflict  with  it.  However,  this  restriction  may  have 
to  be  reevaluated  after  more  work  is  done  using  this  approximation. 

We  will  assume  for  simplicity  that  all  the  sources,  Sk,  are  voltage  sources;  however,  this 
assumption  does  not  mean  a  loss  in  generality,  because  the  two  equivalent  representations  are 
interrelatable.  The  desired  current  /„( jc)  can  be  found  by  summing  up  the  current  response,  /{, 
of  the  cable  n  to  each  source  while  all  the  other  source  voltages  are  set  to  zero, 

=  OHjsjf ;  this  summation  is  then  added  to  the  current  /;  induced  on  cable  n  by  the  incident 


Cable  n  excited  by  the  electric  field  E  and  by  m  —  1  independent  equivalent  sources 
S„(k  =  I ,  nr,  k  *  n),  attached  at  its  ends  (jr„  =  0,  /). 


field,  E,  as  given  by  equation  (1).  Therefore,  the  current  anywhere  along  cable  n  can  be  found 
from 
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where  l0(x)  and  /,( jc)  are  simply  the  current  responses  of  a  transmission  line  excited  by  an  ideal 
voltage  source  of  unit  magnitude  at  each  end,  respectively  (r  =  0  or  /).  Vk  and  Vk  are  ideal 
voltage  source  representations  of  the  source  voltages  Vk  for  x  =  0  and  x  -  I.  The  difference 
between  V*  and  Vk  or  Vk  is  that  V*  is  the  open-circuit  voltage  at  the  end  of  cable  k  and  Vk 
or  Vk  is  the  voltage  at  the  end  of  cable  k  while  loaded  by  the  other  n  -  2  or  m  -  n  -  I  cables. 
Therefore,  Vjf is  related  to  V*!*.,,,-,  by  a  voltage  divider  network  with  the  source 
impedance  Zk  in  series  with  the  parallel  impedance,  ZJ,  of  the  other  n  -  2  source  impedances: 


and 


Similarly,  we  have  for  Vk: 


Vk 


V* 


zt 

zk  +  Zk 


_L  =  "y  _L 

zt  ~  h  z> 


(3) 

(4) 


n  =  v* 


Zi 


Zk  +  Zk 


*«»»+  i.m 


and 


—  =  y  —  I 

Z%  i-»+  |  Zy  I,** 

The  equations  for  l0(x)  and  /,(x)  are1 

/0(x)  =  Kn[Kn  cosh  T „(/,  -  x)  +  ZJ  sinh  rH(ln  -  jc)]/D, 
li(x)  =  -Kn(Kn  cosh  1%*  +  ZJ  sinh  Tnx)/Dn  , 


(5) 

(6) 

(7) 

(8) 


1  R.  F.  Gray.  Nuclear  Kleclromagnetlc  Pulse  Simulation  by  Point-Source  Injection  Techniques  for  Shielded  and 
Unshielded  Penetrations.  Harry  Diamond  Laboratories.  HDL-TR-I7.V7  (December  1975). 


where 


(9) 

(10) 


Dh  =  AC, [(ACJ  +  ZJZF)  sinh  T */„  +  Kn(Zl  +  Zf)  cosh  T „/n ]  .  (11) 

Therefore,  all  the  variables  in  equation  (2)  have  been  defined  by  equations  (3)  through  (11)  except 
for  the  assumed  equivalent  source  terms  Vk  and  Zk  of  the  m  -  1  decoupled  cables.  These  terms 
are  developed  in  the  next  section. 


3.  EQUIVALENT  SOURCE  OF  AN  EMP-EXCITED  CABLE 


The  Thevenin  equivalent  source  terms  can  be  easily  generated  from  equation  (1)  and  the 
geometry  given  in  figure  1.  The  voltage  term  Vk  is  simply  the  open-circuit  voltage  at  the  end  of 
the  Ath  cable  multiplied  by  a  phase  shift  term,  e*""4,  where  ak  is  a  time  delay  which  accounts  for 
the  position  of  the  Ath  cable  relative  to  the  nth  cable.  The  open-circuit  voltage  at  the  end  (jr  = 
0,  /)  of  a  cable  can  be  found  by  multiplying  the  current  i(x)  in  equation  (1)  by  the  appropriate 
terminating  impedance.  Z0  or  Z,.  and  then  finding  the  limit  as  Z0  or  Z,  approaches  infinity: 


V,®.  =  lim  -  Z0/(0) 

Zq-*<X> 

_  E(a>)[(AX  -  jp' Z,)  sinh  f  /  +  (Z,T  -  j/3' K)  cosh  T/  -  Z,(T  + 
(r2  +  /3  '2)[Z,  sinh  H  +  K  cosh  T /] 

Vie  =  lim  Z,i(/) 

Zr*00 


(12) 


E(a>)[(*f  -I-  jfi'Zo)  sinh  Tl  +  (Z0T  +  jfi'K)  cosh  f /  -  Z0(f  -  jp')e10''  ]e~m 
(T2  +  0 ,2)[Z„  sinh  T/  +  K  cosh  f/] 


.  (13) 


It  should  be  pointed  out  that  the  differences  in  equations  (12)  and  (13)  are  due  to  the 
orientation  of  the  cable  with  respect  to  the  incident  electric  field.  If  the  axis  of  the  Ath  cable  is 
parallel  to  the  electric  field  (0*  =  0  deg)  (see  fig.  I)  then  V®„.  =  -  V(K.,  assuming  Z0  =  Zt:  either 
representation  can  be  used  to  find  Vk.  However,  for  dk  ±  0  deg.  the  proper  choice  of  V?K.,  V(„. 
must  be  made  depending  on  the  position  of  the  cable  in  the  electric  field.  If  6k  is  between  0  and 
180  deg.  the  wave  will  excite  the  jr*  =  0  end  of  the  line  first:  therefore.  V?K.  is  used  in  Vk.  If  8k 
is  between  180  and  360  deg,  V'K.  is  required  because  the  wave  hits  the  at*  =  /  end  of  the  line  first. 
The  equation  for  VJh.  has  an  effective  time  shift  for  nonzero  rotation  (180  <  dk  <  360)  because 
of  to  the  phase  constant  /J*.  This  time  shift  must  be  factored  out.  which  is  done  by  multiplying 
V'K,  by  e+w'*'*.  If  0n  is  nonzero,  voltage  terms  will  have  a  time  delay  relative  to  those  at 

x„  =  0;  this  delay  is  accounted  for  in  the  frequency  domain  by  multiplying  the  voltage  spectrum 
by  The  proper  voltage  terms,  V*.  for  either  end  of  the  Mh  cable  and  any  cable  rotation. 

0k,  are  given  in  table  I. 
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TABLE  I 


Voltage  Term,  V*,  for  Equivalent  Source  Sk 


Cable  n 
connection 

0  <  0*  <  180  deg 

180  deg  <  0k  <  360 

.r  =  0 

V“e 

V^e*'*1* 

X  =  / 

The  source  impedance,  Z*,  can  be  found  by  dividing  the  open-circuit  voltage  V(1<.|,_0>,  by 
the  short-circuit  current,  /'sc|x=o.i  which  gives 


0  lim  [Z0/( 0)] 

yn  _  V  m-  _  Zo~»» _ 

*  -i?c  lim  [/(0)] 

Z„-0 


=  K 


K  sinh  n  +  Z,  cosh  T/ 


,  lim  [Z,m 

yl  __  *  <K*  _  0 _ 

*  /'1c  lim  [/'(/)] 


=  a: 


Z,  sinh  T/  +  K  cosh  T/ 

K  sinh  n  +  Z0cosh  T/ 
Z0  sinh  T I  +  K  cosh  T / 


(14) 


(15) 


Equations  (14)  and  (15)  are  simply  the  input  impedances  of  transmission  lines  terminated  in 
impedances  Z,  and  Z0.  respectively.  Obviously,  if  Z0  is  equal  to  Z,  then  Z%  and  Z'k  are  equal. 


This  completes  the  development  of  the  uncoupled  transmission-line  model.  All  the  necessary 
parameters  have  been  defined  and  the  summation  of  currents  on  the  cable  of  interest  (cable  n) 
is  straightforward. 

4.  PENETRATION  THROUGH  AN  IMPERFECT  SHIELD 

If  cable  n  is  actually  a  shielded  cable,  which  is  very  likely  in  tactical  Army  systems,  then  the 
preceding  solution  for  the  current  distribution  on  cable  n  can  be  used  to  calculate  the  internal 
response  of  the  coaxial  cable.  The  response  of  a  coaxial  cable  to  a  single  source  at  only  one  end 
is  given  elsewhere.'  The  solution  for  more  than  one  source  will  simply  result  ir.  a  summation  of 
responses  similar  to  equation  (2).  The  internal  response  of  the  coaxial  cable  due  to  the  current 
coupled  directly  on  cable  n  by  the  incident  field  will  also  have  to  be  included.  The  solution  for 
this  current  (also  given  elsewhere1)  is  part  of  program  FREFLD.  Therefore,  no  additional  devel¬ 
opment  other  than  a  proper  summing  technique  is  needed  for  the  solution  of  the  internal  response 
of  a  coaxial  cable  with  other  cables  attached  to  the  ends  of  its  shield. 

5.  EXPERIMENTAL  VERIFICATION 

A  brief  laboratory  experiment  was  conducted  to  establish  the  validity  of  using  superposition 
techniques  with  transmission  lines.  The  accuracy  of  equation  (1)  has  already  been  addressed.1 
The  overall  validation  of  the  technique,  including  the  assumption  that  the  interaction  between 
cables  is  negligible  for  many  cases,  will  be  investigated  in  future  work. 


1  R.  F.  Gray.  Nuclear  hlectromagnetic  Pulse  Simulation  by  Point-Source  Injection  Techniques  for  Shielded  and 
Unshielded  Penetrations.  Harry  Diamond  Laboratories.  HD1.-TR-I7.17  (December  1975). 
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A  network  of  three  transmission  lines  was  chosen,  as  shown  in  figure  4.  Appropriate  lengths 
of  RG58  A/U  coaxial  cable  were  used  for  the  three  lines.  The  exterior  shields  of  the  cables 
insured  that  the  uncoupled  line  approximation  was  accurate.  The  source  voltages  V,  and  Vt  were 
created  by  placing  small  resistances  (-1  ohm)  in  series  with  the  shields  of  cables  1  and  2  at  the 
junction  with  cable  3.  Two  voltage  sources  were  then  attached  across  these  resistors.  The  source 
impedance  of  the  voltage  sources  was  much  greater  than  I  ohm  so  that  no  loading  of  the  sources 
occurred. 

Three  measurements  were  made  using  an  HP  3042  automatic  network  analyzer  system.  This 
system  consists  of  an  HP  3570A  network  analyzer,  and  an  HP  3330B  synthesizer,  controlled  by 
an  HP  9830  desktop  calculator.  The  calculator  performs  the  data  reduction  and  output  display  as 
well  as  automating  the  measurements.  All  the  measurements  were  made  over  a  frequency  range 
of  100  kHz  to  13  MHz.  First,  the  open-circuit  voltage  V?  at  the  end  of  line  three  was  measured 
relative  to  the  source  voltage  Vt  with  V2  shut  off.  Then  the  voltage  VI  was  measured  with 
respect  to  the  source  voltage  V2  with  V,  shut  off.  Finally,  the  voltage  V3ti  was  measured  with 
respect  to  V,  with  Vt  =  V,. 


v*+ 


Figure  4.  Experimental  test  setup  for  validating  superposition  principle  for  cable  responses. 
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In  order  to  simplify  the  comparison  of  measured  and  calculated  results,  we  assumed  that  the 
transmission  lines  were  lossless.  This  assumption  should  be  relatively  good  except  for  right  at 
the  resonance  points.  The  source  impedances  Z(  and  Z2  are  found  from  equation  (14)  to  be 

.  m,  ,,  -  .  gr  .  2rry/i 


where  K0  =  50  ohms  and  t  =  2x  10*  m/s  for  RG  58  A/U  coaxial  cable.  The  voltage  sources  V, 
and  V2  are  not  needed  since  they  were  equal  in  the  test  and  the  measurements  were  normalized 
to  the  source  voltage  (i.e.,  V3/Vi  or  V3/V2  were  the  measured  values).  The  terminating  imped¬ 
ances  for  cable  3  are  found  from  equations  (9)  and  (10)  to  be 


!/[!/( JKo  tan  0/,)  +  \/(- jK0coi 


The  ideal  voltage  sources  V?  and  V%  at  .v  =  0  are  found  from  equations  (3)  and  (4) 


and  there  are  no  source  terms  at  v  =  I.  The  open-circuit  voltage  VI  at  x  =  /  due  to  either  V ,  or 
V2  is  found  from  equation  (9)  by  multiplying  by  Z(  and  then  taking  the  limit  as  Zt  approaches 
infinity: 

Vfc=  lim  Z,/0(/)  =  Vj|/(ZoCOs  /3/3  +  jK  sin  /3/3)U-i.2  .  (20) 


Therefore,  the  lossless  transmission  line  solutions  for  the  first  two  measurements  using  only 
one  source  are  found  from  equations  (16),  (18),  (19),  and  (20)  to  be 


V2  is  found  from 


The  third  measurement  with  both  sources  exciting  the  system  and  with  V, 
equation  (2)  to  be  the  sum  of  equations  (21)  and  (22). 


The  magnitudes  of  the  experimental  and  calculated  responses  for  the  three  conditions 
(Vj  =  0,  V,  =  0,  Vi  =  Vt)  are  shown  in  figures  5,  6,  and  7.  As  expected,  the  agreement  right 
at  the  resonance  points  is  poor  because  of  the  lossless  line  assumption.  However,  the  comparisons 
of  experimental  and  calculated  results  over  the  remainder  of  the  curves  seem  to  be  more  than 
sufficient  to  validate  the  applicability  of  superposition. 


FREQUENCY  (MHO  FREQUENCY  |MH/| 


Figure  5.  Ratio  of  V?  to  V,  with  Vt  =  0.  Figure  6.  Ratio  of  V?  to  V,  with  V,  =  0. 


FREOu«NCvi»m« 

Figure  7.  Ratio  of  V?,  to  V,  with  V,  =  V,. 

6.  CONCLUSIONS 

An  approximate  solution  has  been  developed  for  the  EMP  response  of  a  coaxial  cable  with 
an  arbitrary  number  of  cables  attached  to  its  ends.  Although  the  effect  of  some  of  the  approxi¬ 
mations  (particularly,  the  noninteraction  of  the  cables)  remains  to  be  determined,  the  basic 
approach  of  superposition  of  transmission-line  responses  has  been  verified  experimentally.  The 
technique  outlined  here  will  be  used  to  develop  a  modified  version  of  the  computer  program 
FREFLD.  Once  this  new  program  is  completed,  the  overall  validity  of  the  technique  will  be 
established. 
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